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METHOD FOR CONSTRUCTING A GEOLOGIC MODEL 

OF A SUBSURFACE RESERVOIR 

[0001] This application claims the benefit of U.S Provisional Patent Application 
5 No. 60/558,287 filed on March 31, 2004. 

FIELD OF THE INVENTION 

'V 

[0002] This invention relates generally to the field of geophysical prospecting. 
More specifically, the invention is a method for constructing a geologic model 
specifying the detailed internal structure and grain size distribution within a 
10 composite sedimentary body based on a seismic image of the body. 

BACKGROUND OF THE INVENTION 

[0003] In the oil and gas industry, seismic prospecting techniques are used to aid 
in the search for and evaluation of subterranean hydrocarbon reservoirs. A seismic 
prospecting operation consists of three stages: data acquisition, data processing, and 
15 data interpretation. The success of the operation depends on satisfactory completion 
of all three stages. 

[0004] In the data acquisition stage, a seismic source is used to generate an 
acoustic signal that propagates into the earth and is at least partially reflected by 
subsurface seismic reflectors. The reflected signals are detected and recorded by an 
20 array of seismic receivers located at or near the surface of the earth, in an overlying 
body of water, or at known depths in boreholes. 

[0005] During the data processing stage, the recorded seismic signals are refined 
and enhanced using a variety of procedures that depend on the nature of the geologic 
structure being investigated and on the characteristics of the raw data. In general, the 
25 purpose of the data processing stage is to produce an image of the subsurface from the 
recorded seismic data for use during the data interpretation stage. 



WO 2005/104003 PCT/US2005/005358 

-2- 



[0006] The purpose of the data interpretation stage is to determine information 
about the subsurface geology of the earth from the processed seismic data. The 
results of the data interpretation stage may be used to determine the general geologic 
structure of a subsurface region, or to locate potential hydrocarbon reservoirs, or to 
5 guide the development of a discovered reservoir. 

[0007] Due to the limited resolution of seismic data, the conclusions which can be 
made after the data interpretation stage are generally limited to broad descriptions of 
the size and overall shape of subsurface reservoirs. The descriptions may for example 
provide an indication of the total volume of hydrocarbons which might be retained in 
10 such reservoirs. However, present technology does not allow the analyst to accurately 
determine the detailed internal structure of the reservoir from seismic data. Even 
when an exploration well has been drilled, present technology does not allow an 
analyst to accurately characterize detailed reservoir structure in locations other than in 
the most immediate region of any such well. 

15 [0008] A commercial hydrocarbon reservoir is generally a "composite" 
sedimentary body, composed of a large number of smaller "fundamental" sand 
bodies, many of which are too small to be separately distinguished in a seismic image. 
Reservoir continuity, the communication (or lack thereof) between adjacent sand 
bodies, is commonly a primary factor controlling hydrocarbon production efficiency 

20 and ultimate recovery. There is a need to predict detailed internal structure of 
subsurface reservoirs using seismic data and without having to drill many exploration 
and delineation wells. Such a capability would facilitate estimation of hydrocarbon 
volume in place and production rates early in the hydrocarbon exploration and 
development process. Accordingly, this invention satisfies that need. 

25 SUMMARY OF THE INVENTION 

[0009] A method for predicting the properties of a composite sedimentary body is 
disclosed. The method comprises (a) determining the outline form of the composite 
body, (b) characterizing the properties of the fundamental bodies in the composite 
body, (c) generating the next fundamental body to be placed in the composite body 
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outline, (d) placing the next fundamental body into the outline of the composite body, 
and (e) repeating steps (c) through (e) until the outline of the composite body is 
substantially full of fundamental bodies. 

[0010] A second embodiment of a method for predicting the properties of a 
5 composite sedimentary body is disclosed. The method comprises (a) determining the 
outline form of the composite body, (b) determining at least part of the outline form of 
at least one fundamental body within the composite body, (c) determining the 
properties of at least one fundamental body within the composite body, (d) 
characterizing the fundamental body properties throughout the composite body, (e) 
10 generating the next fundamental body to be placed in the composite body outline, (f) 
placing the next fundamental body into the outline of the composite body, (g) 
repeating steps (e) through (g) until the outline of the composite body is substantially 
full of fundamental bodies. 

[0011] A third embodiment of a method for predicting the properties of a 
15 composite sedimentary body is disclosed. The method comprises (a) determining the 
outline form of the composite body, (b) measuring the thickness and grain size 
distribution at one point in at least one fundamental body within the composite 
sedimentary body, (c) determining the properties of at least one fundamental body 
within the composite body from the point measurement of thickness and grain size 
20 distribution within the body, (d) characterizing the fundamental body properties 
throughout the composite body, (e) generating the next fundamental body to be placed 
in the composite body outline, (f) placing the next fundamental body into the outline 
of the composite body, (g) repeating steps (e) through (g) until the outline of the 
• composite body is substantially full of fundamental bodies. 

25 BRIEF DESCRIPTION OF THE DRAWINGS 

[0012] Figure 1 depicts the plan view of a fluid flow which is beginning to create 
a fundamental body; 

[0013] Figure 2 depicts the plan view of a fluid flow which has completed 
building one fundamental body and begun building two more; 



WO 2005/104003 



-4- 



PCT/US2005/005358 



[0014] Figure 3 depicts a plan view of a composite sedimentary body containing 
many fundamental bodies; 

[0015] Figure 4 is a flow chart of a first embodiment of the invention; 

[0016] Figure 5 is a flow chart of a second embodiment of the invention; 

5 [001 7] Figure 6 is a flow chart of a third embodiment of the invention. 

DETAILED DESCRIPTION 

[0018] hi the following detailed description, the invention will be described in 
connection with its preferred embodiments. However, to the extent that the following 
description is specific to a particular embodiment or a particular use of the invention, 
10 this is intended to be illustrative only. Accordingly, the invention is not limited to the 
specific embodiments described below, but rather, the invention includes all 
alternatives, modifications, and equivalents falling within the true scope of the 
appended claims. 

[0019] The deposition of clastic sedimentary bodies typically begins with a flow 
of sediment-laden water from a confined channel, such as a river mouth, into an open 
region, such as a basin. Initially such flows expand freely and deposit sediment as the 
flow decelerates. Thereafter, as the deposited sediment grows in height, the deposit 
begins to obstruct the flow. Eventually, the deposit becomes sufficiently large that 
the flow is diverted around it. This terminates construction of the deposit (which we 
call a fundamental body) and results in a new path for the flow to an open region 
beyond or adjacent to it. The deposition process then repeats, and a second 
fundamental body in the system is created. More than one such fundamental 
sedimentary body may be actively built within the system at a time. Overall, the 
process produces a composite sedimentary body consisting of stacks of fundamental 
sedimentary bodies, which is a useful model for the structure of hydrocarbon 
reservoirs. 

[0020] Fundamental bodies are defined as generally lobate sand bodies deposited 
during single flows or multiple flows where the locus of deposition has remained 
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relatively stationary. A fundamental body may also be defined to include other 
elements associated with deposition of the sand, such as the structure and fill 
properties of the channel which supplied the sand to the deposit location, the capacity 
of the flow which built the fundamental body to erode into underlying sediment, and 
5 properties of adjacent fine-grained sediments that may be deposited prior to or 
subsequent to the deposition of the sand. Composite bodies are foimed by aggregates 
! of fundamental bodies and the interaction between elements of the fundamental 
bodies, such as the capacity of the flow which creates one body to erode underlying 
bodies. 

10 [0021] The present invention is a method for constructing a geologic model of a 
subsurface hydrocarbon reservoir using seismic data and, if available, a minimum 
number of well penetrations. The geologic model thus created specifies the grain size 
distribution at all points internal to the reservoir and reflects the fact that the reservoir 
is a composite sedimentary body, made up of smaller fundamental sedimentary 

15 bodies. This geologic model is generated using the outline form of the composite 
body as determined from the seismic data as well as outlines, or partial outlines, of 
some of the fundamental bodies or available well penetrations. 

[0022] In the present invention, the applicants recognized that if the locations and 
properties of the fundamental bodies that build a reservoir could be characterized, 

20 then an accurate geological model of the reservoir could be built. The applicants 
further recognized that due to the large number of fundamental bodies in a 
commercial hydrocarbon reservoir, it is not necessary to exactly specify the location 
and properties of every fundamental body. Rather, any realization of a geologic 
model having the correct statistical distribution of body properties and the correct 

25 trends of change in the distribution with position in the composite body can be used to 
adequately predict hydrocarbon volumes and production rates. 

[0023] Each fundamental body has an associated local inlet, the point at which the 
flow building the body ceased to be confined and became free to expand, as at the end 
of a channel or a local expansion point within a channel. The applicants recognized 
30 that in certain embodiments of the invention, the properties of each fundamental body 
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can be largely characterized by the properties of the depositing flow at the local inlet. 
The grain size trends within a fundamental body as well as the size and shape of the 
body are largely controlled by the flow properties at the inlet, such as the flow 
velocity, flow height, and volume fraction of suspended sediment in each sediment 
5 size class. Additionally, properties of the channel supplying the flow to the body and 
the degree to which the depositing flow erodes sediment within the channel and below 
the new body can be largely determined by the flow properties at the local inlet. 

[0024] The applicants further recognized that the local inlet flow properties for a 
fundamental body may be considered to be a function of a single parameter. The 

10 single parameter may be taken to be the inlet flow velocity, the inlet flow height, the 
inlet Froude number, the inlet width, or any function of those four. When the single 
parameter is determined, the remaining inlet flow properties may be calculated. In 
some embodiments of the invention, the inlet flow velocity is selected as the single 
unknown or free parameter. Assumptions which allow all inlet properties to be 

15 determined from the inlet flow velocity are explained in greater detail in co-pending 
U.S. Provisional Patent Application No. 60/466,655, but are summarized here: (a) that 
the size distribution of grains in the flow is known, (b) that the total volume of 
suspended sediment at the inlet is the maximum carrying capacity of the flow, so that 
the deposition rate is substantially zero at the inlet but becomes non-zero with an 

20 infinitesimal decrease in flow velocity, and (c) that empirical relations giving flow 
height and width as a function of flow velocity for flows in channels apply at the inlet 
location. 

[0025] Figures 1, 2, and 3 depict the process of formation of a composite 
sedimentary body. Figure 1 shows the earliest stage of deposition, where a sediment- 
25 rich flow 10 is emitted from an inlet 12 into an open region 14. As the flow expands 
and decelerates, it loses its capacity to carry sediment, and deposit 16, the first 
fundamental body, is formed beneath the flow. 

[0026] Figure 2 shows a later stage in development of the composite sedimentary 
body initiated in Figure 1. In Figure 2, the first fundamental body, 16 from Figure 1, 
30 has become larger and is labeled 16. Fundamental body 16 is so large that it now 
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diverts the flow 20. Body 16 is no longer building, but the diverted flow is now 
creating two new fundamental bodies, 22 and 24. 

[0027] Figure 3 shows a much later stage in the development of the composite 
sedimentary body begun in Figures 1 and 2. The mapview outline of the composite 
5 sedimentary body 30 and the mapview outlines of two smaller fundamental 
sedimentary bodies 32 within the composite are shown. The inlet of the composite 
body 34 and the local inlets of the fundamental bodies 36 are also indicated. 
Channels, 38, bring the flow from the inlet of the composite body to the local inlets of 
the fundamental bodies. 

10 [0028] A first embodiment of the present invention will now be described, as 
shown in Figure 4. The first step (step 401) is to determine the outline form of the 
composite body. The second step (step 402) is to characterize the properties of the 
fundamental bodies in the composite body. The third step (step 403) is to generate the 
next fundamental body to be placed in the composite body outline. The fourth step 

15 (step 404) is to place the fundamental body into the composite body outline. The fifth 
step (step 405) is to repeat steps 403 thru 405 until the composite body outline is 
substantially full of fundamental bodies. The properties of the fundamental bodies 
include the outline form of the body, the grain size distribution at each point within 
the body, and properties of the flow field that built the body such as the flow velocity, 
20 flow height, flow duration, and suspended sediment volume fractions for each grain 
size range at each point in the depositing flow. The fundamental body properties may 
also include the size, shape, and nature of finer sediment fractions that are deposited 
before or after the main sandy portion of the body, the erosional potential of the flow 
that built the body, and the geometry and fill properties of the channel that conducted 
25 the depositing flow from the composite body inlet to the local inlet of the fundamental 
body. 

[0029] Step 401 : Determine the outline form of the composite body. In one 
embodiment this is accomplished by identifying stratigraphic surfaces of the same 
order that bound, above and below, the composite body. Preferably, the stratigraphic 
30 surfaces extend to the edges of the composite body and close at the edges, though this 
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is not a requirement of the method and partial outlines of composite bodies can be 
used. Methods for identifying stratigraphic surfaces in three-dimensional seismic 
volumes are familiar to persons of ordinary skill in the art, who will also recognize 
other equivalent methods for identifying or inferring the outline form of composite 
5 bodies. Such other methods include interpretation of two-dimensional seismic lines, 
other remote imaging techniques, correlating well logs, and spatially correlated 
outcrop observations. These other methods are also within the scope of this 
invention. 

[0030] Step 402: Characterize the properties of the fundamental bodies in the 
10 composite body. In one embodiment, this step involves characterizing trends in the 
local inlet properties of the fundamental bodies throughout the composite body and 
relating the local inlet properties to all other properties of the fundamental bodies. In 
one embodiment, these trends in local inlet properties can be characterized by 
determining (a) the statistical distribution of local inlet properties at a point within the 
15 composite body and (b) how this statistical distribution changes with local inlet 
position throughout the composite body. The inlet properties characterized in this 
step would be: the flow velocity at the inlet, the flow height at the inlet, volume 
fraction of the flow at the inlet composed of grains within each grain size range, and 
the inlet width. 

20 [0031] In co-pending U.S. Provisional Patent Application No. 60/466,655, 
assumptions are made about relationships between inlet flow properties. In one 
embodiment, flow velocity becomes the only independent variable in defining the 
local inlet properties. Thus, in that embodiment, characterizing the local inlet 
properties throughout the composite body can be accomplished by simply 

25 characterizing the local inlet flow velocities throughout the composite body. 
Typically, the inlet flow velocity trend would be one of exponential decay in local 
inlet velocity with interpreted downstream distance from the composite body inlet and 
the probability distribution would be assumed to be log-normal around the peak 
likelihood specified by the trend. This trend and probability distribution in local inlet 

30 velocity is then fully described by three parameters: the most likely local inlet 
velocity at a point, the characteristic length of the exponential decay in most likely 
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local inlet velocity, and the standard deviation (on a logarithmic scale) of the log- 
normal distribution of local inlet velocities around the most likely value. Methods for 
fitting this three-parameter model to a set of point determinations of local inlet 
velocity are familiar to persons of ordinary skill in the art. One such method is to fit a 
5 straight line to the plot of the natural logarithm of velocity versus distance from the 
composite body inlet. Once the line is fit, the three model parameters can be 
determined as follows: the required velocity at a point can be the velocity at any point 
on the line, the characteristic length of the exponential decay is the negative inverse of 
the slope of the line, and the standard deviation of the logarithm of velocity is the 
10 standard deviation of the measured points around the line. 

[0032] The estimation of parameters of a model describing the variation of inlet 
properties within the system requires actual estimates of inlet properties in at least one 
location within the system. Persons of ordinary skill the in art, with the benefit of the 
disclosures in co-pending U.S. Provisional Patent Applications 60/466,655, 
15 60/459,144 and 60/454,516 will recognize a variety of methods that could be used to 
provide inlet properties for particular sand bodies identified in seismic data or by well 
penetrations. These methods are intended to be within the scope of this invention. 

[0033] Step 403: Generate the next fundamental body to place in the composite 
body outline. In one embodiment, bodies observed in seismic are generated first (by 

20 order of proximity to the system inlet), followed by bodies detected in wells (working 
from the bottom of the well up to the top), followed by enough unobserved bodies to 
substantially fill the composite body outline. To generate a fundamental body, the 
inlet flow velocity must be established. For bodies observed in seismic or in wells, 
the inlet flow velocity is already determined either from the seismic outline of the 

25 body or from the well penetration. For unobserved bodies, the inlet flow velocity is 
chosen at random, based on the probability distribution for inlet velocity at the 
approximate location of the new fundamental body inlet. For the first unobserved 
fundamental body, this approximate location is assumed to be the inlet of the 
composite sedimentary body. For subsequent unobserved fundamental bodies, the 

30 approximate location is assumed to be the inlet location for the previously placed 
fundamental body. Once the inlet velocity for the new fundamental body is 
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established, the other inlet flow properties for the new fundamental body may be 
determined as described in step 402. 

[0034] Once the inlet flow properties for the new fundamental body have been 
established, the fundamental body properties are determined. Typically, the 
5 fundamental bodies are calculated by numerical simulation of the fluid flow which 
deposits them, using the local inlet properties as boundary conditions. A first 
embodiment for making this numerical simulation is described in co-pending U.S. 
Provisional Patent Application number 60/558,287, where time-varying vertically- 
averaged equations for conservation of flow momentum, flow volume, and sediment 

10 volume are solved numerically to determine the evolution of the deposit and its 
internal grain size distribution through time. The simulation progresses until the body 
becomes large enough to substantially divert the flow that forms it. A second 
embodiment for making the numerical simulation is described in co-pending U.S. 
Patent Application number 10/829,600, and is generally preferred over the first 

15 embodiment because it is computationally faster, though less accurate. In this second 
embodiment, vertically averaged equations are also used, but to calculate a steady- 
state flow from the inlet flow properties, assuming flat topography below the flow. 
This flow has an associated deposition rate at each mapview location and is assumed 
to persist in steady state until the body has built to a specified maximum height which 

20 is also a function of the inlet flow properties. Persons of ordinary skill in the art will 
recognize other methods for calculating properties of a fundamental body once the 
inlet flow properties are determined, including other numerical simulation techniques. 
These methods are within the scope of this invention. 

[0035] In many systems, when the flow velocity drops to near zero, very fine- 
25 grained material rains down, creating shale drapes of large lateral extent. These can 
create barriers to hydrocarbon flow in subsurface reservoirs. The thickness and extent 
of the associated shale drape can be calculated for each fundamental body and is 
regarded as a property of the fundamental body. In one embodiment, the approximate 
total volume of shale in the composite body is determined from seismic data using 
30 seismic attribute techniques familiar to persons of ordinary skill in the art. The 
appropriate volume of shale associated with each fundamental body is then calculated 
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by dividing the total volume of shale detected in the composite body by the 
approximate number of fundamental bodies that will be required to fill the composite 
body. The thickness for each shale drape is the volume of shale in the drape divided 
by the lateral area of the drape. Alternatively, a spatially-varying thickness profile for 
5 each shale drape may be assumed which is consistent with the total volume of shale 
and the number of fundamental bodies in the composite body. The size and shape of 
the associated shale drape is recorded as a property of the fundamental body. 

[0036] Underlying sediment can be eroded by the flow fields which deposit 
fundamental bodies. Of particular importance for characterizing reservoir 

10 permeability is the degree to which the shale drapes are penetrated by erosional scour. 
To approximately account for this effect, the flow field and flow duration associated 
with the deposition of each fundamental body is recorded as a property of the body. 
This information will be used later when the fundamental body is placed into the 
composite body outline to determine the extent to which the new fundamental body 

15 will erode underlying sediments. 

[0037] Underlying sediment can also be eroded and new sediment deposited in 
association with the chamiel that brings the sediment from the inlet of the composite 
body to the local inlet of the fundamental body. In one embodiment, the properties of 
the flow in the channel can be assumed to be the same as the properties of the flow at 
20 the inlet of the fundamental body and the channel width may be assumed to equal the 
inlet width. All such channel properties that are derived from the inlet flow properties 
are recorded as a property of the fundamental body. 

[0038] The fundamental body properties, including the associated shale drape, 
erosional flow field, and channel properties, may optionally be pre-calculated for the 
25 full range of inlet flow properties and stored for reference in computer memory to 
reduce the computational time required in application of the inventive method. 

[0039] Step 404: Place the fundamental body into the composite body outline, hi 
one embodiment, this step involves several parts: (a) the location is chosen for the 
fundamental body; (b) the flow field associated with the fundamental body is used to 



WO 2005/104003 



- 12- 



PCT/US2005/005358 



determine the extent to which underlying sediment should be removed; (c) the 
fundamental body is placed into the composite body outline; (d) a channel is cut 
through the composite body model from the fundamental body inlet to the composite 
body inlet, (e) the shale drape associated with the body is placed over the 
5 fundamental body within the composite body outline. These parts are described in 
greater detail below. 

[0040] In one embodiment, the location of the fundamental body is chosen as 
described here. For bodies observed in seismic, the location is already determined. 
For bodies observed in wells, the location of the inlet is constrained by the 
10 requirement that the appropriate thickness contour of the deposit intersect the well. 
The primary (downstream) axis of the fundamental body is constrained to lie on a 
radial line from the system inlet. Subject to these constraints, the local inlet is placed 
as close as possible to the composite body inlet without causing the composite body 
model to exceed the observed thickness of the composite body by more than a 
15 prescribed percentage at any location--the prescribed percentage is typically 25% of 
the fundamental body thickness at the same location. If the location chosen for the 
fundamental body is far from the approximate location used in step 403 to determine 
the inlet properties, step 403 may be repeated using the chosen location as the new 
approximate location for the body. 

[0041] Next, the flow field associated with the fundamental body is used to 
determine the extent to which underlying sediment is removed by erosional scour. A 
variety of empirical relationships between flow properties and erosion have been 
published and are all within the scope of this invention. The preferred embodiment is 
the erosion treatment in co-pending U.S. Provisional Patent Application 60/558,287. 
The appropriate net erosion rate may be calculated at each point below the flow field, 
and sediment removed from the composite body model at each point to a depth equal 
to the erosion rate multipled by the flow duration. 

[0042] Next, the fundamental body is placed into the composite body model. 
Typically, this is accomplished by adding a new layer of cells to the composite body 
30 model, where the thickness and internal grain size distribution in each cell are the 
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thickness and internal grain size distribution of the fundamental body at each 
mapview point. 

[0043] Next, a channel is cut from the inlet of the composite body to the inlet of 
the fundamental body. In one embodiment, the channel is assumed to be of 
5 rectangular aspect, having a depth equal to the inlet flow height and a width equal to 
the inlet width. The sediment removed by this incision is then uniformly mixed and 
redistributed to fill the volume of the incision. Optionally, the channel cut could be 
refilled with sediment having other properties if those properties are known. 
Channels incised based on different assumptions about the cross-sectional shape of 
10 the channel and mapview path it follows from the composite body inlet to the 
fundamental body inlet are also within the scope of this invention. 

[0044] Finally, the shale drape associated with the fundamental body is placed 
into the composite body model. Typically, this is accomplished by adding a new 
layer of cells to the composite body model, where the thickness and internal grain size 
15 distribution in each cell are the thickness and internal grain size distribution of the 
shale drape at each mapview point. 

[0045] Step 405: Repeat steps 403 through 405 until the composite body outline 
is substantially filled with fundamental bodies. 

[0046] A second embodiment of the present invention will now be described as 
20 shown in Figure 5. The first step (step 501) is to determine the outline form of the 
composite body. The second step (step 502) is to identify at least partial outline forms 
of all detectable fundamental bodies within the composite body. The third step (step 
503) is to determine properties of the fundamental bodies. The fourth step (step 504) 
is to characterize the variability of fundamental body properties in the composite 
25 body. The fifth step (step 505) is generate the next fundamental body. The sixth step 
(step 506) is to place the next fundamental body into the outline form of the 
composite body. The seventh step (step 507) is to repeat steps 505-507 until the 
outline form of the composite body is substantially full of fundamental bodies. 

[0047] Step 501: This step is the same as step 401 in Figure 4. 
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[0048] Step 502: Identify at least part of the mapview outline of at least one 
fundamental body within the composite body. This identification is typically made in 
the same three-dimensional seismic volume used to identify the composite body 
outline. The method of picking these smaller units is typically based on contours of 
5 seismic amplitude, but persons of ordinary skill in the art will recognize other 
methods for picking stratigraphic units near the resolution limits of seismic data. 
Automated tools have been described to assist in this process, for example, in co- 
pending U.S. Patent Application No. 10/395,911. Interpreted lines of three- 
dimensional seismic data, other remote imaging techniques, correlating well logs, and 
10 spatially correlated outcrop observations could also be used to identify outlines of 
fundamental bodies within the composite body. These other methods are also within 
the scope of this invention. 

[0049] Step 503: Determine the properties for the observed fundamental bodies. 
One method is specified in co-pending U.S. Provisional Patent Application no. 
15 60/466,655. In this method, the seismic interpreter provides general paleocurrent 
directions sufficient to identify a local inlet point and down-stream end point for each 
identified seismic body. The method then utilizes constraints from fluid mechanics to 
relate a mapview contour of constant deposit thickness to the inlet properties of the 
flow that built the deposit. 

20 [0050] Step 504: Characterize the trend in the properties of the fundamental 
bodies throughout the composite body. In one embodiment, this step involves 
characterizing the inlet flow velocity trend in the fundamental bodies throughout the 
composite body and relating the flow velocity to the other local inlet properties and all 
other properties of the fundamental bodies. As described in step 402, the inlet flow 

25 velocity trend is typically one of exponential decay with interpreted downstream 
distance from the composite body inlet, and the probability distribution around this 
trend could be log-normal around the peak likelihood specified by the trend. This 
trend and probability distribution in local inlet velocity is then fully described by three 
parameters: the local inlet velocity at a point, the characteristic length of the 

30 exponential decay, and the standard deviation (on a logarithmic scale) of the log- 
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normal distribution. These parameters may be estimated from the local inlet 
velocities determined for the identified fundamental bodies in step 503. 

[0051] Step 505: This step is the same as step 403 in Figure 4. 

[0052] Step 506: This step is the same as step 404 in Figure 4. 

5 [0053] Step 507: Repeat steps 505-507 until the composite body outline is 
substantially filled with fundamental bodies. 

[0054] A third embodiment of the present invention will now be described as 
shown in Figure 6. The first step (step 601) is to determine the outline form of the 
composite body. The second step (step 602) is to identify the thickness and grain size 

10 distribution of at least one fundamental body in at least one core that samples 
sedimentary bodies within the system. The third step (step 603) is to determine the 
properties of the sampled fundamental bodies. The fourth step (step 604) is to 
characterize the fundamental body properties in the composite body. The fifth step 
(step 605) is to generate the next fundamental body. The sixth step (step 606) is to 

15. place the fundamental body into the outline form of the composite body. The seventh 
step (step 607) is to repeat steps 605-607 until the outline form of the composite body 
is substantially full of sedimentary bodies. 

[0055] Step 601 : This step is the same as step 401 in Figure 4. 

[0056] Step 602: Measure the thickness and grain size distribution for at least 
20 one fundamental body penetrated by a core sampling the system. Ideally, the core 
should penetrate the composite body itself, but sometimes the lack of nearby cores 
may necessitate using cores that merely penetrate the same depositional system or an 
analogous system. 

[0057] Step 603: Determine the properties of the sampled fundamental bodies. In 
25 one embodiment this involves determining the local inlet conditions associated with 
the sampled fundamental bodies. One method for making this determination is 
described in co-pending U.S. provisional patent application no. 60/459,144, wherein 
the method involves (a) estimating flow properties at the location where the well 
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penetrates each body, (b) extrapolating these flow properties back to the inlet location 
(identified as the point where the flow properties become consistent with flow in a 
channel), and (c) adjusting the initial flow property estimate until the extrapolated 
inlet flow properties are consistent with the extrapolated maximum height of the 
5 body. Applying this method to sand bodies that the well penetrates within the system 
gives inlet properties for these bodies. 

[0058] Step 604: Characterize the fundamental body properties throughout the 
composite body. In one embodiment, this is accomplished by inlet flow velocity 
trend for the fundamental bodies throughout the composite body, and relating the flow 
10 velocity to the other local inlet properties and all other properties of the fundamental 
bodies. As described in step 402, the inlet flow velocity trend is typically one of 
exponential decay with interpreted downstream distance from the system inlet, and 
the probability distribution could be log-normal around the peak likelihood specified 
by the trend. This trend and probability distribution in local inlet velocity is then fully 
15 described by three parameters: the local inlet velocity at a point, the characteristic 
length of the exponential decay, and the standard deviation (on a logarithmic scale) of 
the log-normal distribution. 

[0059] Local inlet properties were calculated in step 603 for bodies in the vicinity 
of the well. If only one well is available, then it may not provide information about 
the trend of those properties throughout the system. One method for determining the 
trend in inlet properties throughout the system uses the relationship between system 
thickness and vertically averaged grain size distribution, as described in co-pending 
U.S. Provisional Patent Application No. 60/454,516. The vertically averaged grain 
size distribution at a point in the system will be approximately equal to the average 
thickness-weighted grain size distribution in the typical sand body at that location. To 
determine the inlet properties of the typical sand body at a given location, the 
following method may be used: (a) determine the vertically averaged grain size 
distribution at that location from the composite body thickness at that location. This 
determination can be made using the thickness and vertically averaged grain size at a 
different location where core is available together with the method of co-pending U.S. 
Provisional Patent Application No. 60/454,516; (b) estimate the typical inlet velocity 



20 



25 



30 



WO 2005/104003 



- 17- 



PCT/US2005/005358 



at that location; (c) calculate the other inlet properties based on the inlet velocity 
estimate; (d) simulate the sand body form and grain size distribution associated with 
the estimated inlet velocity; (e) calculate the average thickness-weighted grain size 
distribution within the simulated sand body; (f) refine the inlet velocity estimate in 
5 step (b) until the average thickness-weighted grain size distribution within the body is 
substantially equal to the vertically averaged grain size distribution at that location as 
predicted in step (a). This method provides the inlet property trends directly from the 
shape of the composite body, without the need for observations of fundamental bodies 
at a second mapview location. 

10 [0060] Step 605: Step 605 is the same as step 403 in Figure 4. 

[0061] Step 606: Step 606 is the same as step 404 in Figure 4. 

[0062] Step 607: repeat steps 605-606 until the composite body outline is 
substantially full with simulated fundamental bodies. 



